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We propose the new method for finding the non-Gaussian tails of probability distribution function 
(PDF) for solutions of a stochastic differential equation, such as convection equation for a passive 
scalar, random driven Navier-Stokes equation etc. Existence of such tails is generally regarded as a 
manifestation of intermittency phenomenon. Our formalism is based on the WKB approximation in 
the functional integral for the conditional probability of large fluctuation. We argue that the main 
contribution to the functional integral is given by a coupled field-force configuration - instanton. As 
an example, we examine the correlation functions of the passive scalar u advected by a large-scale 
velocity field (^-correlated in time. We find the instanton determining the tails of the generating func- 
tional and show that it is different from the instanton that determines the probability distribution 
O . function of high powers of u. We discuss the simplest instantons for the Navier-Stokes equation. 
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^ I I. INTRODUCTION 

> ; 

The intermittency phenomenon (reflected in non-Gaussian, scaling-violating tails of PDF) is believed to be the 
hardest part of the yet to be built theory of turbulence. It is related to rare fluctuations which cannot be treated in 
terms of a perturbation theory. Neither the physical mechanism nor the mathematical properties of such fluctuations 
are known. 

iy~j ' Now, what is the most likely force which can lead to the given rare fluctuation of the field? The main idea of this 
0^ ■ paper is that such force is not random at all. It satisfies the well defined equation, which follows from the WKB 
approximation in the functional integral. Asymptotically, the fluctuations of the force around this most likely one are 
negligible. In this respect, the method is similar to the "optimal fluctuation" method used at treating properties of a 
solid with quenched disorder (see e.g. the book jij) or at treating high-order terms of the perturbation series in the 
quantum field theory 

The problem under consideration is quite general, it can be formulated for any field governed by a nonlinear dynamic 
equation and driven by a random "force" . Generally, the PDF of the field depends both on the statistics of the driven 
force and on the form of the dynamical equation. Here, we are interested in the second dependence so that we assume 
the force to be Gaussian. Because of nonlinearity, the PDF of the field is non-Gaussian even for a Gaussian random 
• 1— 1 , force. Note that strong intermittency appears also for linear problems with "multiplicative noise" , for instance, for a 
' passive scalar advected by a random velocity field. 
We start with the dynamical equation 



d t u + C{u) = <f> , (1.1) 

that controls the evolution of a field u(t, r) under the action of a random "force" cf>(t,r). Here C(u) is a nonlinear 
expression, it can be thought to be local in space. Generally, both the field u and the force (f> have a number of 
components. The Gaussian statistics of the force (p is completely characterized by the pair correlation function 

(^ l! r 1 )^ 2 ,r 2 )) =3(ti-ta,n-r 3 ). (!- 2 ) 



In principle, the relat ions ( LI L^) contain all the information about the statistics of ■ 



The equation (1.1) describes e.g. thermal fluctuations in hydrodynamics where it is reduced to the well known 
Langevin equation |H|. Then (j> is short-correlated in time and in space that is it can be treated as a white noise. 
For some systems, this thermal noise produces remarkable dynamical effects. Some examples are collected in the 
book @. Here, we will be interested in turbulence where 6 is an external "force" correlated on large scales in space. 
Turbulence was first treated in terms of the equation (1.1) by Wyld [0 who formulated the diagram technique as a 



perturbation series with respect to the nonlinear term in the Navier-Stokes equation. The diagram technique cannot 
be applied to our problem since we are interested in non-perturbative effects. Nevertheless we can use the functional 



1 



that generates th e tec hnique since it is a non-perturbative object. Such generating functional was introduced in 
for the equation (1.1) it has the form 



Z(X) 



exp 



dt dr Xu 



T>uDpcxp il 



dt dr Xu 



(1.3) 



where p is an auxiliary field and the effective action is 

X = y dt dr p[d t u + C(u)] 



-- f dtdt'drdr'E(t~t',r~r')pp' . 



(1.4) 



The coefficients of the expansion of Z in A are the correlation functions of u. The auxiliary field p determines response 
functions of the system, for instance, the linear response function (Green function) is G = (up). Note the remarkable 
property 10 



J VuDp exp (il) = 1 



related to causality. That is the reason why the normalization constant is unity in ( |1.3| ). This makes it possible to 
average directly Z over any additional random field if necessary. 

The asymptotics of Z(X) at large A is determined by the saddle-point configuration (usually called classical trajec tory 
or instanton) which should satisfy the following equations obtained by varying the argument of the exponent in (1.3) 
with respect to u and p 



d t u + C(u) = -i j dt' dr' E(t 



t',r-r')p(t',r'). 



dtp - -=-p = X . 
ou 



(1.5) 
(1.6) 



Solu tions of the equations (1.5 L^ ) are generally smooth f unc tions of t and r. Comparing the equations (1.1) and 
(1.5) we conclude that the right hand side of the equation ( |l.5|) just describes a special force configuration necessary 
to produce the instanton. If Ui nst is a solution of ( [l . 5|Jl . 6| ) then asymptotically at large A 

6lnZ(X)/6X = tu inst . (1.7) 



Let us discuss the boundary conditions for the saddle-point equations. The equation (1.5) implies that we should 
fix the value u- m for the field u at the initial time t- m . Contrary, a boundary condition for the field p is implied at the 



remote future since, as follows from (1.6), it propagates in the negative direction in time. Minimization of the action 
generally requires p — > at t — > oo. For the instantons discussed below, the finiteness of the action will also require 
u — > at t — > — oo. 

If one is interested in the simultaneous statistics of u then the function A can be chosen as 



A(t,r) = yS(t)X (v), 



(1.8) 



where y is a number and Ao is an appropriate function of r depending on what spatial correlation function s we are 
going to study. In this case, we should find the solution for p satisfying the rule: p = at t > 0. The system (1.5,1.6) 
is thus to be treated for t < only. That corresponds to the causality principle since only processes occurring in 
the past could influence the value of the simultaneous correlation functions at t = 0. The formal ground for the 
rule follows from the consideration of the problem in the restricted time interval t < to what is possible if A = at 
t > to- Then the minimization of X+ J dtdrXu over the final value u(to) gives p(to) = because of the boundary 
term originating from Jdtdrpdtu. 

One may be interested also in the probability distribution function V(u) for the field u. It can be expressed via the 
generating functional Z(X) via the functional Fourier transform 



V(u) 



J VXZ(X) exp 



dt dr Xu 



(1.9) 



2 



We expect that the behavior of V{u) at large u as well as the behavior of Z (A) for large A i s as sociat ed w ith some 



saddle-point con figurations. Generally, the configurations are not always the same for both (1.3) and (L9). Indeed, 



we see from ( |l.9| ) that the tail of V(u) at large u corresponds to a large value of ShxZ(X)/8X which is related to large 
A only if the tails of both PDF and the generating functional decay faster than exponent - see the example in Sect. 
Otherwise, those tails are determined by different configurations as is demonstrated in Section |n|. 



[II 



The best starting point to develop the instanton formalism is the problem of white-advected passive scalar 9 since 
it allows for a detailed analytical treatment ]l2]-|T^|. It will be shown in Section || that both V(9) and Z{\) have 
exponential tails as it has been established before by Shraiman and Siggia [H| (see also |l4|j ). By using this example, 
we shall explicitly demonstrate that different instanton configurations are responsible for the tails of the generating 
functional Z(X) at large A and of PDF at large 9 respectively. It is instructive to recognize the difference between 
the instantons: We shall show that the instanton that is responsible for large 9 corresponds to a small strain and 
suppressed stretching. Contrary, the instanton that determines the tails of Z corresponds to a large value of strain. 



Section [II presents the first step in studying instantons of the Navier-Stokes equation. Only instantons for the 
two-point generating functional (exp(iA(wi — u 2 )) will be considered. The family of such instantons corresponds to 
the velocity fields with a linear spatial profile at r <C L. Consideration of the instanton perturbations (giving the 
fluctuation contribution into the action) that correspond to spiral creation in the straining field of the instanton will 
be the subject of further publications. 

II. PASSIVE SCALAR ADVECTED BY A LARGE-SCALE VELOCITY FIELD 

Let us show how the general formalism described in the Introduction works for a particular problem: the advection 
of a passive scalar field 8(t, r) by an incompressible turbulent flow in d-dimensional space [[Tl] |l3. The advection is 
governed by the following equation 

(d t +v a V a ~ nA)9 = </>, V a v a =0, (2.1) 

where (f>(t, r) is the external source, v is the advecting velocity and A designates Laplacian, k being the diffusion 
coefficient. Both v(t, r) and (j>(t, r) are random functions of t and r. We regard the statistics of the velocity and of the 
source to be independent. Therefore, all correlation functions of 9 are to be treated as averages over both statistics. 

We assume that the source <p is 5-correlated in time and spatially correlated on a scale L and has a Gaussian 
statistics completely determined by the pair correlation function: 

Mfi,riMt 2 ,r 2 )>=a(ti-f a )x(ria) ■ (2.2) 

Here x( r i2) as a function of the argument r%2 = |rj — r 2 | decays on the scale L. We are interested in the behavior of 
the correlation functions on scales r <C L. Thus, only the constant P 2 = x(0) will enter all the answers. The constant 
P 2 has the physical meaning of the production rate of 9 2 . 

Following Kraichnan jll],[l2| , we consider the case of a Gaussian velocity v delta-correlated in time and containing 
only large-scale space harmonics. Then the velocity statistics is also completely determined by the pair correlation 
function 

(v a (ti,ri)vp(t 2 , r 2 )) = S(ti - t 2 )V a p , 

V a! i = V S a/3 - /C Q/3 (r 12 ), IC a p{0) = . (2.3) 
Here the so-cahed eddy diffusivity is as follows 

IC af3 = D(r 2 S af3 - r a r p ) + D ^ d ~ ^ 8 ap r 2 , (2.4) 



where d is the dimensionality of space and isotropy of the velocity statistics being assumed. The representation 
( pT3 2.4) is valid for the scales less than the velocity infrared cut-off L u , which is supposed to be the largest scale of 



the problem. Then Vq and IC a p in ( |2.3| ) are two first terms of the expansion of the velocity correlation function in 
r/L u so that D ~ Vo/L 2 t . We presume also the inequality dDL 2 ^> k which guarantees the existence of a convective 
interval of scales Td -C r L where correl ation functions of the passive scalar are formed mainly by stretching in the 
velocity field. Here — 2y? n/{D{d — 1)) is the diffusion length 



The statistics of the large-scale velocity field has the remarkable property: It follows from the expressions (2.3,2.4) 
that the correlation function of the strain field a a p = V pv a is r-independent: 

(<rap(h)<7nu(te)) =D[(d+ l)5 afi 5p v ~ 5^5$^ - 8^5^] 5{t\ - t 2 ) , (2.5) 
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That means that the strain field a a p can be treated as a random function of time t only. Just that property enables 
one to find in detail statistical properties of the field 9 Jl^Jl^]. To exploit the property, it is convenient to pass into 
the comoving reference frame that is to the frame moving with the velocity of a Lagrangian particle of the fluid. That 
means that we pass to the new space variable r — g(t) where g(t) is the Lagrangian trajectory of the particle H[|. 
We will take the particle positioned at the origin at time t — 0, then 



g(t) = J dTv(T,g{ T )) 



(2.6) 



After the transformation r — > r — g(t), the equation (2.1) acquires the form 

{d t + [v a (t,r) - v a (t,0)]V a - kA} = <j> 



(2.7) 



It is seen from (2~3|,2.4) that the statistics of v a (t, r) — v a (t, 0) coincides with the statistics of u a prp. That means 
that the generating functional corresponding to (|2.7|) can be written as 



Z{\) = J V9VpVcrexp \ + i j dtdv\9 



(2.8) 



where a a p is a function of time satisfying a aa — 0, its PDF is exp(— J-). The effective action I and the functional T 
in (§J) are 



il = i J dtdr (pd t + pa a prp\7 a 9 + n\7pS79) 
-i J dtdr-i dr 2 Pix(ri2)p2 , 

dt [(d + l)<T a pcr a p + o- a/3 a/3 a ] 



T = 



2d(d + 2)D 



(2.9) 
(2.10) 



Note that there is the difference between (1.3) and (2.8) which is in the presence of an additional random field a a p. 



A. Uniaxial Instanton 



Here we examine the saddle-point contribution to the generating functional Z (A). The equ ations d escribing the 
saddle points are extremum conditions for iT + i J dtdr X9 — T. Starting from the expressions (2.9 2. 1C ) we find 



d t 9 + a a(3 r p V a 9 - nV 2 9 = -i J dr' X (\r - r'\)p(t, r') , 
d t p + (TaprpVap + nV 2 p = A , 

cr a p(t) =iD dr ((d + l)r /3 V a 9 - r a Vp9 - S a pr 7 V 7 9^p , 



(2.11) 
(2.12) 
(2.13) 



where p = p{t, r), 9 — 9(t,r). If to take into account only the saddle-point contribution described by ( 2.11 , 2.12 , 2.13 ) 
then 



Z{\) = ^exp (i J dtdr X9^j \ cx cxp(— JF oxtr ) , 
where ^"oxtr is the saddle-point value of T — il — i J dtdr \9. One get from (2.9, [2.10 2.12) 

^"extr = ^ J dtdr 1 dr 2 pix{ri2)P2 



(2.14) 



2d(d + 2)D 



dt [(d + l)a a pa a p + a a pap a ] 



(2.15) 



In t he f ollowing, we will be interested in simultaneous correlation functions of 9 so that we tak e the field A in the 
form (1.8) and solve the equations for negative time t < 0. Let us stress that for the function (1.8) the term X9 is not 
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influenced by the t rans formation r — > r — g(t) because of g(Q) = 0. Note that the system of equations (2.11,2.12,2. 13 ) 
with the function (1.8) is invariant under the transformation 



a — > Xo , p — » Ap , t — > X 1 t , y ^ Xy , k — > Xk , ^-" e; 
where X is an arbitrary factor. It leads to the conclusion that 

•^extr = yf{y/K) , 



XT c: 



(2.16) 



(2.17) 



with the function / to be determined. 

We will treat nearly single-point statistics. That means that the space support of the function Ao in fll.8Q is taken 
to be much smaller than the pumping length L. From the other hand, we would like to avoid bulky formulas related 
to the account of diffusion. Therefore, the size of the support is believed to be much larger than the diffusion length 
Td- We thus come to 



X(t,v)=y6(t)S A (v), 



(2.18) 



where <5a(i") is a function with the characteristic size A 1 satisfying L ^> A 1 ^> ra and normalized: J dr^A(r) = 1. 
For example, we can take 



Sa(t) 



A d 

Zdj2 



exp(-A 2 r 2 ) . 



We thus examine the following object 



where 



Za = (exp(iy9 A )) , 



9 A = J dr* A (r)0(t = O,r). 



(2.19) 



(2.20) 



(2.21) 



Having in mind the in equality A 1 3> r^, we omit in the following the diffusive terms in the equations. The extremum 
conditions (2.11,2.12) are then as follows: 



'xp'- 



d t + (J a (3rp\7 a 9 — -i J dr 
d t p + VaprpVaP = yS(t)5 A (r) , 



(2.22) 
(2.23) 



It is natural to seek a solution of (2.22,2.23,2.13) in the uniaxial form what means that a a p is a diagonal matrix 
with the components 



diag a = (— s, s/(d — 1), . . .) . 
As it was suggested in , it is useful to pass to the new fields 

0(t, r) = 6(t, e«x, e±r±) , p(t, r) = p(t, e«x, e±r±) . 



(2.24) 



(2.25) 



where x is the coordinate along the marked direction, is the component of the radius-vector r perpendicular to 
the direction and 



ei|(t') = exp 



dts(t) 



d-l -1 

e i = e„ 



Now, the equations (2.22,2.23) can be rewritten as 



d t 6 = -i J dr' X (R(t))p(t,r') , 

d t p = yS(t)5 A (r) -> p = -yS A (r) , t < , 



(2.26) 



(2.27) 
(2.28) 



where we presented an obvious solution for p satisfying p — at t > 0. The quantity R in (2.27) is 
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R 2 = e 2 (x - x') 2 + e 2 ± (r ± - r' ± ) 2 . 

Note that 

d t R = -s{xd x + x'd' x )R + ^( ri Vi + r' ± V ± )R . 
For the considered uniaxial geometry, the relation ( p. 13 ) gives 

s = -iD J drp[(d- l)xd x - r±V±]0. 

Using now fl2.27U2.28Q we find 

d t s = -Dy 2 J dvdv 1 S A (r)S A (v') [(d - l)xd x - r± V J X (J2) . 

By virtue of (2.30) and the symmetry properties of the integrand in ( 2. 32] ), we obtain 

(d-l)Dy 2 



(2.29) 
(2.30) 

(2.31) 

(2.32) 

(2.33) 
if t -> -oo 

(2.34) 

One can demo nstrat e that the main contribution to Z A is determined by the saddle point with s > 0. The n, in 
accordance with ( 2.26 ), en increases with increasing \t\. That means that the characteristic value of R in ( 2.34 ) can 
be estimated as R ~ A. _1 e||. At small \t\ where is not very large, x(-R) in ( 2.34 ) can be substituted by P2 = x(0) 
and we find that s ~ si where 



s9*s = 



ft / drdr'£ A (r)5 A (r')x(#). 



The equati on ha s an obvious first integral which can be established if to take into account that s 
[otherwise fl2.15 ) is infinite]: 

s 2 = (d-l)Dy 2 ( dvdv'8 A {v)5 A (v') X {R) . 



Sl = yy /(d-l)P 2 D 



(2.35) 



That leads to e m « exp(si|i|), which is correct if R < L what means \t\ < s 1 1 ln(LA). In the opposite limit 
|t| ^> s^ 1 ln(LA), the value of s tends to zero. 

The above analysis shows that the main contribution to .Foxtr (2.15) is associated with the region of integration 
|i| < s^ 1 ln(LA) and the first term in fl2.15| ) can be written as {l/Z)y 2 P 2 s^ 1 ln(LA). Substituting ( p.24| ) with ( |2.35| ) 
into the second term of (2.15), we find 



To 



P2V 2 



(d-l)D 



ln(LA) . 



(2.36) 



Note that the expression is in agreement with ( 2.17 ) since we considered the case r^A <C 1 where the answer should 
be K-independent. It is also possible to restore Z A (y) in the limit A — > 00 that is for the single-point object. F or thi s 
we should recognize that generally ^"extr is a function of the dimensionless parameter Ar^ and use the property ( 2.17 ). 
Then in the limit r^A ^> 1 where A-dependence should disappear we obtain 



3~ 2: 



' P2V 2 
(d-l)D 



\n{L/r d ) + i In 



P 



2 2 



D 



y\d-i) 



(2.37) 



Note the nontrivial dependence of this single-point object on y. It is the consequence of the time-dependence of the 
effective diffusion cut-off which could be seen at the direct solution with an explicit account of the diffusion. 

Above, we considered the simplest case of the uniaxial strain matrix a a p. It is not very difficult to generalize the 
scheme for the case where principal axes of <j a p are fixed (that is do not depend on time) . The answer shows that it 
is the uniaxial solution that gives the minimum value of .Fextr and therefore only this contribution should be taken 
into account. 

As long as w e are interesting in the tail of the generating function Z A {y) at large y, the instanton contribution 
(2.36) or (2.37) gives the correct answer. However, it is not enough to consider that contribution to obtain the tails 
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of the PDF because the respective tail of Z(y) is exponential. Indeed, we shall see below that the tails of V(9) are 
determined by the contributions at moderate y. We thus face the problem of finding 2(y) at arbitrary y. Fortunately, 
the tails of V(9) are also determined by the instanton contribution yet this instanton is different from the above 
uniaxial solution, which represents the situation where stretching occurs along one marked direction. It is obvious 
that if the direction slowly varies in time the value of the effective action will not be essentially influenced. The role 
of such soft fluctuations is expected to be negl igible if the characteristic time s^ 1 ln(LA) of the stretching is small 
enough. We thus conclude, taking into account ( 2.35| ), that the expression ( 2.36| ) is correct at large y. At moderate y, 



the fluctuations of the stretching direction should be taken into account, it is the topic of the next subsection. There, 
we shall explicitly integrate over soft mode and obtain different equations for the instanton. 

B. Isotropic Instanton 

Here, we are going to take into account the fluctuations of the stretching direction which were neglected in the 
preceding subsection. For that purpose, it is useful to introduce the variable measuring the stretching rate along 
the current stretching direction (the direction of the maximal Lyapunov exponent) determined by th e stra in field 



Pap = V/3i> Q . For this aim, it is useful to perform the transformation of the fields 9, p generalizing (2.25) for an 
arbitrary a a p [Q . Namely, let us pass to the fields 

8{t, r) = 9(t, M afi rp) , p(t, r) = p(t, M aP r ) , (2.38) 
with d x d matrix M a p controlled by the equation 

d t M = bM , M(t = 0) = i , (2.39) 

with a formal solution 

t 

M 



Texp N dtf&(tffj . (2.40) 
o 

The symbol T designates the anti-chronological ordering for negative t. Note that detM = 1 due to the incompress- 



ibility condition tr<r = V Q w Q = 0. Performing the substitution in (2.9) and passing to the new space variable Mr 



(the Jacobian of the transformation is equal to unity due to detM = 1), one obtains 

%L = i J dt drpd t - i J dt dr x dr 2 p x x{R)p2 , (2.41) 

where 

R a = M aP (r w - r 2 p) . (2.42) 

We see that only R is er-dependent in ( [2.4l| ) and, moreover, only it's absolute value R enters the effective action. 
Just that value is a measure of the stretching irrespective of the directions of the current main axes of the matrix a. 
The statistics of R can be established starting from the PDF exp(— T), see e.g. jl5|. The answer is that, for negative 
times, R can be written as 



R(t) 



exp(y dt'((t'))\r 1 -v 2 \, (2.43) 



with the random variable £ having PDF exp(— J-q) with 



The generating functional (2.20) is thus rewritten as 

Z A = [ V9VpVQ exp (iyO A + il - T c ) , (2.45) 
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where 9\ is defined by ( [2.21 ). 

We have performed the exact transformation of the statistical weight. Let us stress that the situation described by 
( [2.45 ) is isotropic from the beginning whereas the solution found in the preceding subsection was aniso tropic. Now, 
we formulate extremum conditions describing a saddle point for the argument of the exponent in (2.45): 



dtp = y5(t)5 x (r) , 

d t 0{t,T X ) = -i I dr 2X [R{t)}p(t,r 2 ), 



«t>) = d Jtp±D- d -^±D J dtdr^p^R. 



(2.46) 
(2.47) 

(2.48) 



(2.45 2.46) that in the saddle-point approximation Z\ cx exp(— .Tvjxtr) where 



The equat ion (2.46) has the same form as (2.28) and has consequently the same solution p = — y8A(r). It follows from 

(2.49) 



If 1 

•^cxtr = ^ J dtdr l dr 2PlP2X(R) + 2D{d- 1) 



It follows from fl2.43| ) that d t R = -£R. Using that, we can find the first integral of the equation (E 



d 2 {d~\f 



D 2 + (d - l)Dy 2 J dri dv 2 8 K ^i)^2)X{R) ■ 



(2.50) 



The c onstant here is established using the property £ — * (d — l)dD/2 at t — > — oo f ollow ing from 9 — ► at t — * — oo 
[( |2.4£ ) is infinite otherwise]. The characteristic value of R in the right-hand side of ( 2.5CQ can be estimated as 



R(f) 



A _1 exp / dt((t) . 



(2.51) 



If R <C L then the integral in the right-hand side of (2.5C) is approximately equal to P 2 , if R 3> L then the integral 
is negligible. That means that there are two different time intervals. At large \t\, it is C — (d — l)dD/2 and at small 
\t\ it is C — Ci where 



2 (d-l) 2 d 2 n2 



D z + {d- l)Dy z P 2 



(2.52) 



The boundary between the regions is at \t\ - Ci"" 1 MLh). The main contribution to JFcxtr fl2.49| ) is associated with 
the region |t| <h = Cf 1 ln(LA). The first term in (|2.49[ ) can be substituted by y 2 P 2 t\/2 and the second one can be 
substituted by (2D(d- l)) _1 (Ci - (d- l)d£»/2) 2 i x . Using ( ^52| ) we find 




P 2 y 2 



(d-l)D 2 



ln(LA) . 



(2.53) 



Comparing the expression ( 2.53| ) with ( |2.36| ), we conclude that fluct uatio ns of the stretching direction can be neglected 
if y 2 Dd z P^\ Let us stress that at y 2 ~ Dd^P^ 1 the value of ( ^.53 ) is much larger than unity. That means that 
violation of (2.36) is not associated with destructing saddle-point regime, it is rather related to an incorrect calculation 
of soft fluctuations in the saddle-point regime. Note also that the role of fluctuations increases with increasing the 
space dimensionality d. 



C. Probability Distribution Function 

T he sc heme proposed in the preceding subsections can be applied also to calculating PDF Va($) of the quantity 
6 \ ( 2.21 ) . Let us start from the average 



9|") 



Jv8VpVaexp(il-F a + 2n\n8A). (2.54) 
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Thus we sec that the saddle-point contribution to (9 A n ) is determined by extrema of iT — T a + 2n In 9\ which coincide 



with (pl]j2.12| , ^13|) if to substitute 



2ni 
~0~A~ 



(2.55) 



Then an attempt to find the analog of the uniaxial instanton fails. The formal reason for this is in additional i in 
(|H|). The physical reason is that the uniaxial instanton is an adequate object for the statistics of fast processes 
whereas (9 A n ) is determined by slow processes. 

To find a solution, we should pass to the isotropic instanton. That means that we should per form the same 
transf orma tion of the fields as in the preceding subsection what leads to the saddle-point equations ( 2.46 , 2.47 , 2.48 ) 
with ( [2.55 ). The equations has a solution of the same type as considered above with 



d 2 {d~lf 



D 2 -(d- 1)DP : 



An 2 
2 — 

A 



(2.56) 



The value of 9a in (2.56) is the parameter which can be found from the equation analogous to (2.47) which now reads 

2n 



dtO 



Oa 



dr 2 x(#)<5A(r 2 ) 



(2.57) 



As previously, the integral in the right-hand side of ( 2.57 ) for ri = is equal to P 2 if \t\ m (£A-)/d and is negligible 
otherwise. We thus come to the conclusion that 



\t = 0,r = 0) ~ 2nPi 



ln(LA) 



Substituting the relation into ( |2.56 ) we find the equation on d leading to 

Ci = (d-l)D\- 



ln(LA) V In 2 (LA) 



(P 

T 



(2.58) 



(2.59) 



We see that Ci decrea ses wi th inc reasin g n and consequently the characteristic time ln(LA)/£i increases with increasing 
n. Substituting now ( [2~59l) into ( ^58|) one obtains 



9 2 

'An 



8nP 2 ln(LA) j n n 2 d 2 

d 2 D{d - 1) i MLK) + V In 2 (LA) + T 



(2.60) 



It is not very difficult to recognize that the main contribution to the saddle-point value of iT — + 2n In 9a is 
determined by the last term. That means that 



oc exp(-J r cxtr ) oc 9 



An ' 



(2.61) 



with 9 An from fl2.60| ). 

The same result can be deduced by t he al ternative method. Namely, starting from ( 2.53 ) we can calculate the tail 
of the PDF "Pa(i9) for the quantity 9a (2.21). The function Va^) is the Fourier transform of Za{u)- 



Pa(i?) = J dy exp(-iyd)Z A (y) oc J dy exp(-iy& - T^) . 
Substituting here ( [2.49] ) and calculating the integral over y by the saddle-point method |ll| we find 



7> A (tf)ocexp^ -ln(LA) I 1- Wl 



d- 1 



Pi 



-D- 



& 2 



\n\LK) 



(2.62) 



(2.63) 



what is in agreement with |T3|,|l5j,yT| . Formally, the expression (2.63) is valid at $ — > oo but really it covers the whole 
region of d because the PDF is Gaussian at small $ |1J] . The distant tails of the PDF are exponential as has been 
established first by S hraiman and Siggia |l3| . Note that the value of the Lyapunov exponent (4 corresponding to the 
saddle point in (2.62) is 
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d ( d - 1 $ 2 \ 1 

C oxtr = -(d -1)D[1+ . (2.64) 

S 2 V ' \ P 2 In 2 (LA) J V ' 



We see that the value decreases with increasing ■& whereas the value Ci ( 2.52 ) increases with increasing y. Note also 
that the value of y corresponding to the extremum point is 



., {d-l)d 2 D tf 2 



VcxtI 4P 2 tf 2 + P 2 /{{d-l)D) ' ( 2 ' 65 ^ 

That means th at |?/ 2 Yt ,.| < d 3 D/ P 2 and conseq uentl y the extremum point lies beyond the applicability region of the 



approximation ( 2.36| ). This is the reason why (2.36) does not admit to restore Va($) 



Now, we can calculate {0\ n ) starting from the definition 

dd$ 2n V K {$). (2.66) 



9f) 



This integral can be calcu late d usin g again the saddle-point method. The result coincides, of course, with ( 2.61 ). We 



thus conclude that ( |2.49 ) or (2.63) cover both cases of slow and of fast processes. That means that the account of 



fluctuations in the direction of stretching (performed in the preceding subsection) give us the tool for finding tails of 
both PDF and of the generating functional. 

D. Discussion 

We considered the statistics of the passive scalar advected by the random velocity field in the framework of the 
instanton formalism. The consideration was very instructive since it revealed some nontrivial peculiarities of the 
formalism. First of all, we see that the direct solution of the saddle-point equations gives us the answer which 
satisfactory describes the tail of the generating functional Z(X) but cannot serve to restore the tail of PDF V^ff). 
The physical reason for this lies in difference of processes forming the tails: The tail by Z is related to the fast 
processes with the characteristic time decreasing as A increases while the tail by V is related to slow processes with 



characteristic time increasing as d increases. The conclusion can be directly extracted from (2.35) and (2.64). For 
the slow processes, the fluctuations of the stretching direction are relevant which do not destroy the saddle-point 
(instanton) regime but renormalize the naive answer. For that particular problem, the fluctuations can be explicitly 
taken into account after the special transformations of the fields. Although the trick cannot be widely generalized 
it shows the direction of improving naive answers. In the general case, we expect that the direct solution of the 
saddle-point equations will produce nonsymmetric instantons with a degeneracy parameter (like the direction of the 
marked axis in the considered case). Then, there exists the "Goldstone" mode related to slightly nonhomogeneous 
variations of the parameter. Such mode is soft since it only weakly disturbs the naive instantons. Therefore, the 
fluctuations related to the soft mode are relevant and should be explicitly taken into account. 

III. THE SIMPLEST INSTANTON OF AN INCOMPRESSIBLE VELOCITY FIELD 

Here, we describe the first step in considering a much more complicated problem of finding the tails of the PDF 
for velocity field in three-dimensional incompressible turbulence. We consider two-point statistics and show that an 
instanton with a linear spatial profile naturally appears as a basic flow. 



The effective action (1.4) for the Navier-Stokes equation can be written as follows 



T = 



J dt dr (p a d t v a + p a vp\7f3V a - vp a \7 2 v a + p a \7 a P + QV q d q ) 

J dtdt'drdr'E(t-t',r-r')p a p' a , (3.1) 



The additional independent fields P and Q play the role of Lagrange multipliers enforcing the incompressibility 
conditions V q Uq, = and the analogous condition V a p a = for the response field p a . The field P has the meaning 
of pressure (divided by the mass density p). The origin of the terms with the fields P, Q in the effective action is 
related to the continuity equation dtp + V a {pv a ) = 0, which should be incorporated into the effective action, Q is 
just the auxiliary (response) field corresponding to the equation. At the condition that all velocities are much smaller 
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than the sound velocity, it is possible to ne glec t the time derivative in dtp + V Q (pw Q ) = and variations of the mass 
density what leads to the term QV a v a in (3T). While variations of the mass density can be neglected variations in 
the pressure are relevant. Therefore, it is natural to pass from the integration over the mass density to the integration 
over the pressure as it is implied in (3.1). 

We are going to examine the generating functional for the velocity 



Z{\) 



exp 



dt dr Av 



VpVvVPVQ exp il 



dt dr Av 



(3.2) 



The extremum conditions for the argument of the exponent in (3.2) determining the Navier-Stokes instanton read 

d t v a (r) + vp(t,r)VpV a (t,r)-vV 2 v a (t,r)+V a P(t,r) = -i J dt' J exp(ikr) S(i - t', k)p a (t', k) , (3.3) 

d t p a (t, r) - pp(t, r)V a vp(t, r) + vp{t, r)Vpp a (t, r) + v^ 2 Va (t, r) + V a Q(t, r) = \ a (t, r) , (3.4) 

where 3(k) and p a (k) are Fourier transforms of 3(r) and p Q (r) respectively. In ( |3~3|3~4 ) the conditions V ' a v a = 0, 
VqPq, = are also implied which originate from varying over the fields P and Q. Then the values of the fields P and 
Q can also be found from the conditions. It gives the relations 



V 2 P = -V a (vpVpv a ). 

V 2 Q = V a (jppV a vp - vpS7pp a ) . 



(3.5) 
(3.6) 



Note that similar idea of the instanton formalism has been discussed also by Giles JTt| yet his approach is based on 
uncontrollable approximations. 

In the following we consider the simultaneous correlation functions of the velocity differences ([v(0,p/2) — 
v(0, — p/2)] 2 ") where p is the separation between the points. The functional generating such functions is extracted 
from Z(X) if one gets 



X a = yn a 6(t)[S(v - p/2) - 5(r + p/2)] , 



(3.7) 



where n is a unit vector. As was explained in Introduction the presence of such term in right-hand side of (3.4) means 
that we should solve the problem at negative times t with the final condition 



Pa (0, r) = -y(5 a p - V a VpV- 2 )np[S(r - p/2) - <5(r + p/2)] 



(3.8) 



We assume that th e pumping correlation function S is delta-correlated i n tim e: 3(t, r) = S(t)x(r). Then the system 
of equations (3.3-3.6) is invariant under the transformation analogous to (2.17) 



X~H. 



Av, P^X 2 P, 



X 3 Q, 



Xv , A — > XX , 



X 3 p, 



(3.9) 



where X is an arbitrary factor. For the function (3.8) the transformation ( [3.9] ) means y — > X 2 y. The extremum value 
•^cxtr of the argument of the exponent in (|3^) transforms as .Fextr — * A" 3 J r 0Xtr at ([}!]). That leads to the conclusion 
that 



Z(y) oc exp(-J r oxtr ) J^xtr = y 3/2 f(y/v 2 ) ■ 



(3.10) 



with the function / to be determined. We expect that in the limit y — > oo a ^-dependence in the function disappears. 
Then we conclude .Fextr oc y 3 ^ 2 . 

The characteristic wave vector kg in the correlation function x(k) of the pumping force is of the order of the inverse 
pumping length L. Then examining the region r <C L one can expand the exponent exp(ikr) in (3.3) into the series 



over kr. The first term of the expansion produces the zero contribution to the right-hand side of (3.3) because of 
the structure of the field p determined by the condition (3.8): The condition means that at t = p(r) = — p(— r), 
the property is reproduced by the equations, so that p(k) — p(— k) at any tim e t. Thus the leading term of the 
expansion of the right-hand side is linear in r. That means that the equation ( |3.3| ) admits as a solution in the region 
|r| -C L a linear profile 



cr a p(t)rp , a aa = . 



(3.11) 
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Then we obtain from (3.3) 



d t CT al + OafiOfi-t — -5 



d d k 



(2tt) 



^fc 7 p a (k)x(k) . 



Here we substituted the expression for the pressure 



P 



1 



(3.12) 



(3.13) 



which provides for the condition V Q t) a = 0. Note that P is defined up to a harmonic function, the expression ( 3.13| ) 
is chosen because of its isotropy. 

For the linear velocity profile the equation ([O]) can be rewritten in Fourier representation as 

d 

dtPa - crpaPp - apikfi-^-p a - vk 2 p a + ik a Q = (3-14) 



which should be solved with the condition following from (3. 



p a (t = 0, k) = 2iy S a/ 3 



k a k^ 



np sin(kp/2) . 



(3.15) 



The characteristic wave vector k in J3.12j ) is of the order of L~ l . Thus we can expand sin(kp/2) in kp and keep 
only the first nonvanishing term of the expansion cx kp. As was d iscus sed in the Introduction, the response field 
p(r, t) propagates backwards in time, starting with the initial value (3.15) at t = 0. We shall see that for a long time 
(determined by a small viscosity) the field p(t, k) at k ~ L^ 1 has t he sam e structure cx kp. 

There is a general family of the flows with linear profiles - see Sect. [II C below. We start by considering the simplest 
case. We assume below that the point separation p is directed along the same vector n a as the measured velocity 
components: p a = n a p. Then the problem possesses the axial symmetry which allows us to look for the following 
uniaxial strain matrix 



a al3 = s (6 af3 - dn a np) 



The same symmetry admits the anzatz 



p a (t,k) = 5 a 



k a k^ 



iynp<j>(t, z)knp , 



(3.16) 



(3.17) 



correct in the limit of small k. In (|3.17 ) z = kn jk and the function <f>(t, z) to be found has the initial (final) value 
4>(t = 0, z) = 1. Substituting Q3.15|j3.17D into (pU^ ) we find Q = and 



l (d t - vk 2 )(j> + dz(l - z 2 )d z (f> + 2[(d - 1) - dz 2 ](f) = . 



Substituting the expression (3.17) into the right-hand side of (3.12) we find 

d d k 

f fc 7 p a (k)x(k) = G{dn a n 1 - S ai ) , 



(2tt)' 



(3.18) 
(3.19) 



where 



G = iypC / dzz 2 {\ - z 2 )( d - 1)/2 cf)(t,z) 



(3.20) 



and 



C 



Sd- 



(2Tt) d {d- 1) 



dk k d+1 X (k) . 



(3.21) 



Here Sd is the area of the unit sphere in (i-dimensional space Sd — 2it d / 2 /T(d/2). The const ant C can be estimated 
as C ~ £/ L 2 where £ = (v a dtv a ) is the energy dissipation rate. Substituting now j3.16tl3.19D into (|I|) we find 

d t s = (d-2)s 2 -G. (3.22) 

Our next problem is to find G as a functional of s, to close this set of equations. We have to solve the equation ( 3.18| ) 
for (p. Here, viscous and inviscid cases are slightly different. We start from considering an inviscid Euler equation, 
then we will account for the viscosity. 
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A. Instanton of the Euler equation 



Neglecting viscosity in (3. IS) we get a general solution 

■ z 2 h -2d 



h 2 z~ 2 F 



1-z 2 



h(t)=expfj s(t')dt' 



The initial condition <fi(0, z) = 1 fixes the function F: 



,2~2d 



(*) = 



1 - z 2 + z 2 h~ 2d 



We get the system of equations 



= {d-2)s 2 -G{h), h = sh. 



(3.23) 



(3.24) 



(3.25) 



This system for the variable q — h 2 becomes the usual potential problem q = —U'{q) with the potential 



For d = 3, 



G(h) 



U(q) = -(d-2) I dqqG(q£*) . 

2 1 \n(2h 2d - 1) 



iypCh 2 
h 2d -l 



1 - h' 2d 3 



h 2d -l 



The relevant solution, which vanishes at t — — oo, corresponds to the zero energy in this potential [s(t) oc —1/t as 
t -> -oo]. Therefore, q 2 = 2U(0) - 2U(q) and 



g t=0 = y/2[U(0)-U(l)] = dy/Syp/L. 



(3.26) 



Then, the strain at the moment t = becomes a a p = q(6 a p — dn a np) / 'q{2 — d). In accordance with (1.7) the log- 
arithmic derivative of Z-functional is related to the average initial value of the velocity difference Z'(y)/Z(y) = 
n a {v a {p,0) — v a (—p,0)). In the leading WKB approximation at large y this average can be replaced by the contri- 
bution from the instanton solution: 



(lnZ)'(y) = 2pn a np<T a f3 = — ^- — — - = C 2 \/ '£yp 3 L 2 . 



q(d-2) 



Finally, we obtain the surprisingly simple result 



Z{y) oc exp E{ypfL- 2 



(3.27) 



(3.28) 



with the dimensionless constants C\, Ci and C3 to be calculated. This result is in agreement with the general form 
( |3.10 ), it contains also p-dependence. 



B. Account of viscosity 

When the viscous terms are kept, the solution is modified as follows. With the same assumption k p ~ p/L -c 1, 
we can still look for the unifo rm strain solution. The viscosity will drop from the velocity equation, but not from the 
response field equation (3.14). The extra term vk 2 p can be compensated by extra time dependent exponential 



/;,,(/. k) = j 8 af} - ) iynp<f>(t, z)knpexp [uR(t, z)fc 21 



(3.29) 



The balance of k° and k 1 terms in the equation is the same as before. The balance of k 2 terms gives the equation 



f) /? 

R=l + sLR = 1 - dsz(l - z 2 )-^; + 2s(l - dz 2 )R 



(3.30) 



13 



with the boundary condition R(0) — 0. The substitution R(z,t) = A(t) + B(t)z 2 reduces the PDE to two ODEs 



,4 = 1 + 2As 

B = 2s(B - Bd- Ad) . 

The solution is expressed via s(t); at t — > — oo it grows linearly: Rm t. The influence of the viscosity on our solution 
is weak, it smears the peaks at p and manifests itself when vR ~ L 2 i.e. at t ~ L 2 jv. That time should be much 



larger than the time of instanton formation ^/L 2 jEyp. Our asymptotic expression ( 3.2S| ) is insensitive to viscosity if 



C. Instanton family 



Considering more general strain does not change basic conclusions of this section. Let us describe, for instance, a 
general three-dimensional symmetric flow of the type considered in Jl8| . In the cylindrical coordinates with z-axis 
along p, the velocity vector field at r <C L is given by 

u = (u r ,ug,u z ) = (—ar/2,ujr/2,az) . (3.31) 

Here vorticity has only z-component u/{t) which is a function of time as well as strain cr(t). The pressure is now of 
the form 

P = —gr 2 — e [r 2 sin 9 cos 9 + rz (sin 9 + cos 9)] . 

Particular details of the solution depend on the relation between g and e. The diagonal elements (proportional to 
g) are determined locally from Poisson equation AP = —div(u ■ Vu). Note that the off-diagonal pressure elements 
are generally determined by the global structure of the flow. In our case, the value of e is determined by the distant 
asymptotics u — ► 0, P — >const at r — > oo and matching conditions at r ~ L, which depends on the particular choice 
of the pumping \- The global description of the flows for the whole instanton family is still ahead of us. As far as the 
functional dependence of the respective Z(y, p) is concerned, it is the same for the whole family and does not depend 
of the large-scale behavior of the pumping. Considering, for instance, the case 5 = [opposite to the diagonal case 
( [3-13 ) considered before], we get u> = y/3s and the system of equations similar to ( 3.25 ) 



..2 



S 



G'(h), h = sh 



with another yet qualitatively similar function G' . For the variable q(t), related to s(t) by q = —sq, the Newton 
equation appears with a potential energy U that allow for a single solution (zero energy separatrix) vanishing as 
s(t) o c l/t at t — > — oo. The basic result \nZ(y,p) oc (yp) 3 ^ 2 ^/£ /L is valid for the whole family in agreement with 



(3.10) 



D. Discussion 



The particular instanton found has the scaling 

5u(p) = u(p/2) - u{-p/2) oc p , (3.32) 

it would give the asymptotics of the right tail of the PDF T'(8u,p) oc exp[— (5u/p) 3 ]. It is unclear at the moment if 
there are flows where such an asymptotics takes place; most probably, this simplest instanton does not realize the main 
extremum of the action. Note that the similar instanton with the linear profile is found for the Burgers problem |T^ ] 
where it indeed determines the right tail (Su > 0) of velocity PDF due to sawtooth waves. The general analysis of the 
whole family of instanton solutions for the two-point velocity statistics at the framework of the Navier-Stokes equation 
will be published elsewhere. Also, the crucial problem of the contribution to the action from the fluctuations against 
the instanton background will be considered. It is clear that, in the straining flow of the instanton, any vorticity 
perturbation produces a spiral with the accumulation point at the velocity null. The scaling of the perturbation 
contribution is different fro m (|3.32| ); for instance, it will give Kolmogorov's 5/3-law for the pair correlation function 
as in the Lundgren example |19||. The analysis of the instanton fluctuations will be the subject of further publications. 
Note that the instanton formalism provides a natural (and long-expected) tool for incorporating numerous results on 
particular solutions of the Navier-Stokes equations into the statistical theory of turbulence. 
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